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Abstract - Droplet formation of suspensions is present in many industrial and technological processes such 
as coating and food engineering. Whilst the finite time singularity of the minimum neck diameter in capillary 
break-up of simple liquids can be described by well known self-similarity solutions, the pinching of non- 
Brownian suspension depends in a complex way on the particle dynamics in the thinning thread. Here we 
focus on the very dilute regime where the filament contains only isolated beads to identify the physical mech¬ 
anisms leading to the pronounced acceleration of the filament thinning observed. This accelerated regime is 
characterized by an asymmetric shape of the filament with an enhanced curvature that depends on the size and 
the spatial distribution of the particles within the capillary thread. 


Droplets are present in our daily life and their formation 
has been investigated scientifically since many decades or even 
centuries. Yet, the physical mechanisms of capillary break-up 
of even simple fluids has been only understood within the last 
20 years dEl. The final stages of break up are governed by 
self similar laws that depend only on the fluid parameters such 
as surface tension, viscosity and density ED. The situation is 
less clear for complex fluids such as polymer solutions Ed or 
suspensions 0 - 11 ]. Suspensions can be found in many indus¬ 
trial applications, such as painting, food processing or cosmet¬ 
ics and a better understanding of the capillary break up process 
will be important for a better control of e.g. the dosage. 

Non-Brownian suspensions have recently attracted interest 
as their detachment dynamics are function of the individual 
and collective particle dynamics in the thread and cannot be 
described using an effective fluid approach only. Several dis¬ 
tinct regimes have been identified during suspension detach¬ 
ment CMS and in particular, an acceleration of the detach¬ 
ment compared to the interstitial fluid has been observed during 
the very late stages of droplet detachment fl6l . Numerical sim¬ 
ulations m have confirmed these experimental findings and 
have also indicated the exitance of an accelerated regime where 
the thinning rate is faster than the thinning rate of the pure sol¬ 
vent. While this acceleration has been experimentally quanti¬ 
fied as a function of the number of particles in the thread m 
the exact thinning dynamics in the accelerated regime were so 
far outside the experimental resolution. Here we present mea¬ 
surements on dilute suspensions where only very few particles 


or even no particles at all are present in the thinning capillary 
bridge. At these small concentrations the viscosity of the sus¬ 
pensions can be considered to remain equal to the viscosity of 
the suspending fluid and the change in thinning dynamics can 
unambiguously be attributed to the presence of single beads in 
the thread. High speed imaging in combination with a ’’super 
resolution technique” G3 are used to characterize the break up 
dynamics (Fig. |T]). 

A capillary bridge of sufficient length is prone to the 
Rayleigh Plateau instability. This is also true for a droplet 
hanging from a faucet EHH, and the capillary neck first thins 
exponentially in time with a characteristic time constant that 
depends on the diameter of the neck and the fluid parame¬ 
ters. Close to the break up, the dynamics are described by 
self similar laws, and depending on the fluid parameters, one 
first observes a regime that is governed by viscosity (the Stokes 
regime) and then by viscosity and inertia (the Navier-Stokes 
regime) d. Recently, it was shown that in the transition be¬ 
tween these two regimes a short lasting pure inertial regime 
might be identified, and of course, for other fluid parameters, 
different sequences can be observed (20|. Here we show that 
the presence of even single particles in or near the capillary 
neck strongly perturbs the thinning dynamics and leads to a new 
regime with an accelerated thinning of the filament. Very close 
to pinch-off, the thinning dynamics eventually return to the 
self-similar Navier-Stokes regime, when the filament is fully 
depleted of beads. We quantify for the first time the transition 
points and the velocity of this accelerated regime and quantita- 
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Fig. 1: High speed images of capillary break-up of a 3 wt% suspension 
of 80 pm particles (top sequence) and silicon oil (bottom sequence). 
Frames are taken at the same minimum neck diameters, i.e. at times 
before pinch-off t = x p — i = 17.53, 12.33, 5.19, 3.90, and 1.29 ms for 
the suspension and t - 19.48 , 15.58, 7.79, 5.84, and 1.94 ms for the 
oil before pinch-off. The scale bar is 0.2 mm. 

tively show that they depend on the size and the spatial distri¬ 
butions of particles in the capillary neck. 

Our capillary break-up experiments are performed in a 
home-made capillary extensional break-up rheometer (CaBER) 
device. A well defined quantity of the sample liquid is placed 
with a pipette between two steel discs with a diameter of 2 mm. 
The lower disc is fixed and the upper disc can be drawn apart 
very slowly with a linear motor (P01-23x80, Linmot, Spreit- 
enbach, Switzerland). To measure the minimal neck diameter 
h(t) as a function of time, a LED lamp is used to produce a 
shadowgraph image of the capillary bridge. Light that passes 
through the center of the filament is not diffracted and pro¬ 
duces a bright line in the middle of the filament, while the light 
that interacts with the boundaries of the filament gets diffracted 
so they appear dark. In principle, the intensity profile of the 
diffracted light can be evaluated to achieve spatial resolutions 
below lOOnm ED- The filament is filmed by a 10 bit high 
speed camera (X-Stream XS-5, IDT, Tallahassee, USA) at up 
to 8.2 kHz frame rate depending on the spatial resolution. The 
camera has 1280 x 1024 pixels with a size of 12 x 12 jum. The 
camera is equipped with a tube lens and microscope objectives 
(Nikon) with magnifications ranging from 2 to 10 fold. 
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Fig. 2: (Color online) The minimum neck diameter h(t) of the thin¬ 
ning capillary bridge of the silicon oil (circles) and the suspensions 
(triangles). The data shown for the silicon oil is identical on all three 
panels. From left to right: 80 pm((|) = 3%), 40 pm((|) =1%) and 
20 pm((|) = 3%). The red lines have a fixed slope of 0.015ms -1 and in¬ 
dicate the Stokes regime. The blue lines have been fitted with a linear 
law and indicate the accelerated regime (see also Fig. [6]). The thinning 
dynamics close to pinch-off for the area indicated by the square on 
panel (a) are shown in Fig. fusing experiments with a higher resolu¬ 
tion. 

Our samples have been described before m. In brief, 
spherical polystyrene beads Dynoseeds from Microbeads with 
grain diameters d = 20,40 and 80 pm at volume fractions (|) = 
Vg/V o = 1 to 3%, with V g the volume of grains and Vo the total 
volume, are dispersed in a silicon oil (Shin Etsu SE KF-6011) 
closely matching their density. The viscosity of the pure oil 
is T|o = 0.18 Pa.s and its surface tension a = 21 db 1 mN/m at 
T = 21° C. The volume fractions were chosen such that in most 
experimental runs only single beads would be present in the 
capillary bridge. 

At these low volume fractions, the viscosity of the suspen¬ 
sions is well described by Einsteins formula r\ = T|o(l + 2.5<|>). 
The maximum difference between the viscosity of the suspen¬ 
sion and the interstitial fluid is between 2.5% or 7% and this 
small change in viscosity has negligible effects on the thinning 
dynamics of the viscous thread he 

A typical series of images close to pinch-off for a sample 
with 80 pm beads at cf> = 3% is shown in Fig. QJtop). One 
observes that single beads remain in the thinning thread, lead¬ 
ing to perturbations of the free surface and a larger curvature. 
The distribution of the beads becomes inhomogeneous as the 
thread thins further and at some moment the thinning dynamics 
become localized between two beads and finally the capillary 
bridge breaks between the two beads that are furthest apart. In 
the following we will show that the exact thinning dynamics 
of this regime depend on the particle size and the individual 
distribution of beads in the thread. 

If we define the time t = T p — x as the difference of the actual 
time x from the moment of pinch-off x p and plot the minimum 
diameter of the capillary bridge h(t) as a function of time for 
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Fig. 3: (Color online) The minimum neck diameter h(t) of the thinning 
capillary bridge of the silicon oil (circles) and the 80 pm suspension 
(triangles) corresponding to the results within the black square in Fig. 
13 (a) obtained from experiments with a higher resolution. Lines with 
fixed slopes of 0.007ms -1 (black) and 0.015ms -1 (red) illustrate the 
transition from the Stokes to the Navier-Stokes regime for the oil. For 
the suspensions, a linear fit yields a slope of 0.02ms -1 (blue) for times 
t > 2ms. 


the pure silicon oil (Fig. [2 black circles) we recognize the well 
known different regimes of capillary thinning. First, for times 
t » 10 ms before pinch-off, the dynamics are still governed 
by the exponential growth of the primary Rayleigh Plateau in¬ 
stability (data not fully shown). Approaching pinch-off, the 
system will then first follow the viscosity dominated Stokes (S) 
self similar law m 
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With higher temporal and spatial resolution (black circles in 
Fig. 0, a further transition to the Navier-Stokes (NS) self sim¬ 
ilar regime with 
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is observed at t^s ~ 2ms corresponding to a minimal diame¬ 
ter of h^s ~ 14 pm. 

For the suspensions (Fig. [2 black triangles) the slopes of 
the Stokes regime are identical to the slopes for the pure oil, 
confirming our hypothesis that at the small volume fractions 
we work with, the suspension viscosity can be considered to 
be identical to the viscosity of the interstitial fluid. However, 
the capillary thread does not evolve directly from the Stokes to 
the Navier-Stokes regime regime for the suspensions, but we 
observe an intermediate accelerated regime with a higher slope 
compared to the Stokes regime (Fig. [2. The average slope of 
this accelerated regime for all our experimental runs is 0.022 d= 
0.002 m/s. No significant dependence on the particle size is 
observed within our experimental resolution, but the slope of a 
given experiment depends on the distribution of particles within 
the filament as we will show below. 

At even later stages of the detachment the Navier-Stokes 
regime, with a slope that is identical to the pure oil, is recovered 
again (Fig. [2). 
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Fig. 4: (Color online) Shapes of the thinning capillary filaments. The z 
axis represents the vertical direction with gravity being oriented from 
left to right, (a) For the oil (b) for the suspension with 80 pm beads 
(<|> = 3%) (b). The shapes corrospondend to times t = 18.18, 14.28, 
12.33, 10.39, 7.79, 5.84, 4.54, and 1.94 ms.(c) One minus the corre¬ 
lation coefficient (the residual sum of squares) of a parabolic fit of the 
shapes, indicating the symmetric shape for the pure oil (rhombi) and 
asymmetric for the suspension (squares). 


The Stokes and the accelerated regime are not only distin¬ 
guishable by their different scalings of the minimum neck di¬ 
ameters but also by their symmetries along the z-axis. The 
shape of the thread in the Stokes regime is z-symmetric around 
the minimum diameter ID, as is shown for the pure oil in Fig. 
0] a. In contrast, at the same minimum filament diameter, the 
filament shapes of the suspension have not only a much larger 
curvature but they evolve also from a symmetric shape in the 
Stokes regime to a highly asymmetric shape in the accelerated 
regime (Fig. 0]b). This is quantitatively shown by represent¬ 
ing the correlation coefficient (the residual sum of squares) of 
a quadratic fit to the filament shape (Fig. 0]c). The asymmetric 
profile of the thread also proves that the accelerated regime is 
not a transient inertial regime that has recently been observed 
in Newtonian fluids l20l and where the profile is symmetric. 
However, it is the specific shape of the filament of the suspen¬ 
sion with its large curvature that leads to the accelerated regime 
and we will describe in the following how the size and distri¬ 
bution of the particles affect this regime. 

The minimal neck diameters at which the crossover to the ac- 
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Fig. 5: (Color online) a) Histogram of transition diameters h a for the 
suspensions from Stokes to the accelerated regime. The dashed line 
indicates the transition diameters h^s at the transition from the Stokes 
to the Navier-Stokes regime for the pure oil. b) Histogram of the dis¬ 
tribution of transition diameters from the accelerated regime to the 
Navier Stokes regime h^s- The dashed line indicates the transition 
diameters h^s at the transition from the Stokes to the Navier-Stokes 
regime for the pure oil. 


celerated regime and subsequently to the Navier-Stokes regime 
occur are function of the particle size, as can already be seen 
from Fig. [2] This indicates that the transition is caused by 
the perturbation of the thread due to the presence of beads (see 
Fig. ID- To quantify this further, we have measured the minimal 
neck diameters at these transition points for a total of 236 ex¬ 
periments for the three different particle diameters. Histograms 
of the minimal diameters at the transition to the accelerated 
regime are shown in figure \5\ a. It is clearly observed that for 
smaller bead sizes the transition takes place at smaller neck di¬ 
ameters. For the smallest bead size the transition to the ac¬ 
celerated regime sometimes overlaps with the transition to the 
Navier-Stokes regime for the pure oil. The increase of the tran¬ 
sition diameter with particle size is in agreement with previous 
observations by Bonnoit et al. ifTHl who obtained average tran¬ 
sition diameters comparable to those reported here. They also 
showed that the dependence of the average transition diameter 
on the particle size does not follow a simple scaling law and 
seems to be slightly less than linear. Here we show that there 
is a large variation in transition diameter for a given bead size, 
as can be seen from the histograms of Fig. \5\ We will show 
below that this large variation can be explained by the varying 
distribution of individual particles in the viscous thread close 
to pinch-off. Fig. \5\ b shows the minimal neck diameters at 
the crossover from the accelerated to the Navier-Stokes regime. 
Once more, the transition diameters h^s become smaller with 
decreasing particle diameter and they are always smaller than 
the transition diameter h^s between the Stokes and the Navier- 
Stokes regime for the pure oil (Fig. 0b). As before, the transi¬ 
tion diameters h^s show variation for a given particle size due 



Fig. 6: (Color online) a) Highspeed images for different experimen¬ 
tal runs at « 7 ms before pinch-off (accelarated regime). The arrows 
indicate the distance D between the two beads where final break-up 
will occur. The slope of the minimum diameter vs. time curve h(t) 
are from left to right: 0.025, 0.024, 0.023 and 0.022 mm/ms. b) The 
slope of the accelerated regime for the 80 pm(6 = 3%) suspension as 
a function of the distance between the two beads that are farthest apart 
in the filament at t = 1ms before pinch-off (cp. Fig. 1). The line is a 
guide for the eye only. 


to the variation in particle distributions in the thread. 

We will now quantify the acceleration caused by the par¬ 
ticles. To do so we have measured the slope in the acceler¬ 
ated regime and have linked it to the largest distance between 
two particles in the viscous thread (Fig. [6]). In the acceler¬ 
ated regime the filament thinning is localized between these 
two particles and the viscous thread will break there. We chose 
to measure this distance at 7 ms before the pinch-off where 
the thinning dynamics are well into the accelerated regime. We 
have performed this analysis for the largest particle size leading 
to the highest precision when determining particle positions in 
the thread. If we look at the snapshots at t = 7 ms we find that 
the filament will always break between the two particles that 
are most apart and that the slope of the h(t) curve is a function 
of the distance D of these two particles (Fig. ia>. The smaller 
this distance is, the larger is the curvature and the dynamics will 
be fastest. Interestingly, we find that the accelerated regime is 
observed even if the beads are only at the ends of the capil¬ 
lary thread, in agreement with observations by van Deen et al. 
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m. The thinning dynamics in the accelerated regime depend 
thus strongly on the microscopic distribution of the particles. 
The varying particle distribution is most likely also at the ori¬ 
gin of the large variations in the values of the minimum neck 
diameter at the crossover to and from the accelerated regime. 
Note that a theoretical approach by Hameed and Morris ll22l 
predicts a slowing down of the thinning of a viscous thread in 
the presence of particles, but in their analysis the particles were 
held fixed contrary to the experimental situation, probably at 
the origin of the observed differences. 

In conclusion we have shown that the final stages of pinch- 
off in a very dilute suspension are accelerated compared to the 
pure suspending fluid, despite the fact that the bulk viscosity of 
the suspensions is nearly unchanged by the presence of parti¬ 
cles at these small concentrations. We show that the filament 
thinning dynamics for the suspensions evolve through a self¬ 
similar regime identical to the one of the intersticial fluid be¬ 
fore entering an accelerated regime. At even later stages of 
the detachment process the suspensions recover the self-similar 
solution of the interstitial fluid and the thinning dynamics be¬ 
come identical to the interstitial fluid again. The acceleration is 
caused by the perturbaion of the thread by the particles in the 
viscous thread. They locally introduce a curvature that is too 
large to accommodate the self similar solution of the Stokes 
regime. This leads to the strong acceleration of the thinning 
dynamics until the self similar solution is recovered and the 
thinning dynamics become identical to the pure fluid again. 

The local distribution of isolated particles captured in the 
thread determines the filament diameter at which the crossover 
to the accelerated regime takes place. As the exact distribution 
varies from one experimental run to another a large variation 
is observed for the crossover diameters. The average value is 
function of the particle diameter, but no simple scaling between 
the crossover filament diameter and the particle size has been 
found. Similar observations are made for the filament diameter 
at which the accelerated regime crosses over to the self-similar 
solution again. The slope of the accelerated regime is also a 
direct consequence of the distribution of particles in the thread 
and we have shown that it is function of the distance between 
the last two remaining particles in the filament. 

We have thus shown that the presence of isolated particles in 
the viscous thread strongly affect the thinning dynamics close 
to final pinch-off. The exact evolution of the shape of the vis¬ 
cous thread in time is function of the particle size and distribu¬ 
tion in the thread. 

* * * 
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